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O ' ABSTRACT 

u 



It is discovered in our previous work that different observational systematics, 
e.g., errors of antenna pointing directions, asynchronous between the attitude and 
science data, can generate pseudo-dipole signal in full-sky maps of the cosmic mi- 
crowave background (CMB) anisotropy published by The Wilkinson Microwave 
Anisotropy Probe (WMAP) team. Now the antenna sidelobe response to the 
Doppler signal is found to be able to produce similar effect as well. In this 
CN) ' work, independent to the sources, we uniformly model the pseudo-dipole signal 

in ' 



and remove it from published WMAP7 CMB maps by model fitting. The result 



CN | demonstrates that most of the released WMAP CMB quadrupole is artificial. 

in 

q ' Subject headings: cosmic microwave background -- cosmology: observations 



X ' 1. PSEUDO-DIPOLE SIGNAL 

a: 

The WMAP mission makes measurements of temperature with two antennas, A and 
B, and record in time-order the difference between the two antenna temperatures, Ta — Tb, 
which is called the time-order data (TOD). The observed CMB signal is contaminated by 
Doppler effect induced by the joint motion of the solar system and the spacecraft. The 
aroused dipole difference signal can be calculated by 

d= — v-(n A -n B ), (1) 

c 

where T = 2.725 K is the CMB monopole, c is the speed of light, v is the joi nt velocity, n A 



and n B are the unit direction vectors of the antenna A and B respectively fjHinshaw et al. 
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20091 ). In this paper the vectors v, n A , and n B are defined in the spacecraft coordinate 
system. The errors of line-of-sight (LOS) vector, An A and An B , will produce a pseudo- 
dipole difference signal 

d* = ^v(An A -An B ) 

= — v-An. (2) 

c 

The dipoles for each observations have to be removed from the raw data before map- 
making because their intensities are roughly 10 to 20 times greater than those of the CMB 
anisotropies. A small error of antenna direction will produce an error in predicted dipole 
intensity and then cause a pseudo-dipole signal in the resulting CMB map more noticeably 
through the Doppler dipole subtraction process^ . For example, a LOS error of ~ 7', just 
about a half-pixel in the WMAP resolution, can consequently cause the dipole signal to be 
deviated by 10-20 /iK, which can not be ignored compared to the very weak CMB signal. An 
asynchronous between the attitude and differential data can also produce the pseudo-dipole 
signal. The WMAP mission uses two separate clocks for the attitude data and science data 
respectively. Therefore, if there is a small constant timing error, there will be a constant 
direction difference between the "observed pixel" a nd the true pixel. This has the same effect 



as a constant LOS error in spacecraft coordinates (JLiu. Xiong fc Lill2010l ). 



2. SIDELOBE PICKUP 

Another possible source of pseudo-dipole signal in released WMAP maps is the sidelobe 
signal contamination. Like all radio telescopes, the WMAP antennas have both main beam 
response and sidelo be response. The WMAP antenna sidelobe response was described by 



Barnes et al.l (120031 ) and the corresponding data file is publicly available!. The data files are 
fits format full sky maps in spacecraft coordinates in which the sidelobe responses are given 
in normalized gain G, where the normalization rule is that the summation of all gains for 
one antenna (including the main beam) equals to N, the number of pixels in the map. Thus 
for each differential observation, the recorded difference signal is X^=o i^t ~~ ^f )^V^- Let 



*A true and pure time-order dipole signal produces only a dipole component in the final CMB map, but 
a time-order pseudo-dipole signal can produce more components, especially quadrupolc, on the final CMB 
map. See Fig. 1. 

2 The data file of WMAP sidelobe response can be found at 



http://lambda.gsfc.nasa.gov/product/map/dr4/farsidelobeinfo.cfm 
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k denotes pixels corresponding to a certain sidelobe, then the overall sidelobe response is 
12k(^k ~ Gk)Tk/N. The overall sidelobe response should contribute a pseudo component 
to the differential Doppler signal 

d* - T °y y ( g ^~ g fc B ) n k m 

a sidelobe ~ v /__, yy • l°J 

it 

The spacecraft LOS vectors and n k are all constant vectors in spacecraft coordinates. Since 
G are normalized gains, it can also be treated as constant. Therefore we have two constant 
vectors: An^ = J2k ^k n k/^ an< ^ ^ n B = S/c ^k n k/^- Since the sidelobe gain is much 
smaller than the main lobe gain, An^ and An^ are two small constant vectors, and 

d* siMobe = ^v.(An;-An;) 

= ^v-An*. (4) 

c 

By comparing Eq. 2 and Eq. 4, we can see that the entire sidelobe Doppler pickup can be 
exactly described by a small constant deviation to the differential spacecraft LOS vector. 
Similarly, the effect of the sidelobe response uncertainty on the differential Doppler dipole is 
also equivalent to introducing a small differential LOS error An*. It is important to notice 
that, even if the LOS is absolutely accurate, this pseudo-dipole signal from the sidelobe 
response uncertainty can still exist. 

The amplitude of the equivalent LOS error induced by sidelobe Doppler pickup can be 
estimated by arcsin(|An*J/|n A |), which are ~ 50' — 75' for different bands. The WMAP 
team believes that the Galactic sidelobe pickup is neglectable except for the K-band, the 
reason is largely tha t the Galactic emission is much weaker in other bands than in the K-band 



([Barnes et al.ll2003l ). The case is significantly different for the Doppler signals because they 
have almost the same amplitude in all bands. Moreover, the Doppler signal is as strong as 
several mK all over the sky; however, the Galactic emission is strong only for the low Galactic 
latitude and the strength decreases rapidly for higher latitude. Thus the overall power of the 
Doppler signal is about ~ 10 — 40 times higher than the overall Galactic foreground emission 
power, and the effect of sidelobe Doppler signal pickup is at least 10 times stronger than the 
sidelobe Galactic pickup. 

The WMAP antenna sidelobe gain patterns are estimated by ground-based measure- 
ments and in-flight lunar measurements. For the K-band, the WMAP first-year in-flight 
gain measurements are 60% systematically brighter than ground-based measurements. The 
WMAP team scaled up the ground-based measurements by 30% and scaled down the lu- 
nar results by 30% to yield a best-guess sidelobe gain map with an overall calibration un- 
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Fig. 1. — Induced map from the pseudo-dipole signal produced by a differential direction error An, for 
WMAP7 year-1 Ql-band, in mK and in Galactic coordinates. From left to right: St x with An = (0.01,0,0) 
, Sty with An = (0,0.01,0), 5t z with An = (0,0,0.01). Results for other years and other bands are similar 
to this. 



certainty ~ 30% (jBarnes et al.ll2003l ). thus the amplitude of the equivalent LOS error in- 
duced by the sidelobe Doppler pickup uncertainty could be up to ~ 15' — 22' (estimated 
by arcsin(|An A |/|n A |) x 30%), which is much more than enough to produce considerable 
pseudo-dipole contamination. From the WMAP releases we find that the WMAP7 sidelobe 
file is exactly the same with WMAP1 for the K-band, and for other bands the sidelobe files 
of a single survey year are just slightly changed from year-1 to year-7, indicating that their 
uncertainties have not been essentially improved and can not be ignored. 



MODELING AND REMOVING PSEUDO-DIPOLE SIGNAL 



The approach we use to remove the pseudo-dipole signal is same to the WMAP fore- 
ground removal. Different effects arising foregroun d are separately mode led with template 
maps: the 94 GHz dust map t d for dust emissio n (JFinkbeiner et al.lll999l ). the full-sky Ha 
map t H for free-free emission (jFinkbeinerl 120031). and the synchro t ron template t. 9 d erived 



from K, Ka bands for synchrotron emission (IBennett et al. 



2003 



foreground emission is removed by the WMAP team (IBennett et al. 
with model fitting to get the clean temperature 



2003 



Goldet al.ll2009l\. The 



Gold et al.ll2009h 



Mclean £ ' \Cd *d ~t~ ^h h ' ^s^s) j 

where the coefficients Cd, c H , cs are determined by weighted fitting to make t 



(5) 



2 
clean 



mm. 



The pseudo-dipole signal in WMAP temperature maps can be removed with the same 
approach for foreground contamination. The overall effect of the three systematics mentioned 
in §1 and §2, LOS error, time drift, and sidelobe Doppler pickup, can be generally described 
by an equivalent deviation An of the differential antenna direction. The deviation induced 
by the pseudo-dipole signal in a released WMAP map can be easily modeled. For each 
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measured temperature difference in the TOD used to produce the map, we substitute the 
pseudo-dipole signal calculated by Eq. 2 for an assumed An to produce a new TOD|j . 
The temperature map produced by map makings from the new TOD can be used as a 
model of the pseudo-dipole signal in the released map. In calculation we use the spacecraft 
coordinate system (X,Y,Z), where the X axis is parallel to plane of radiators, — Z is the 
anti-sun direction of the spin axis, and Y is perpendicular to both. The WMAP spacecraft 
scans the sky with a hybrid motion mode consists of rotation and precessing. In spacecraft 
coordinates, the LOS unit vectors of its two antennas are close to (x, y, z) = (0, 0.94, —0.33) 
and (x, y, z) = (0, —0.94, —0.33), and the spacecraft rotation is around the Z-axis. Suppose 
the overall LOS error An from all such effects is made up of three small vectors (5,0,0), 
(0,5,0) and (0,0,5) with 5 = 0.01, each alone on the final map induces its own full-sky 
distribution of deviation St x , 5t y and 5t z , respectively. The induced deviation upon the 
released WMAP7 year-1 Ql-band map, St x , St y , and St z are shown in Fig. 1. Results for 
other years and other bands are similar. 

It is easy to see from Fig. 1 that St y and St z are highly correlated. In order to avoid 
degeneracy issue we only use 5t x and St y in model fitting. The clean full sky temperature 
map 

tdean = tf — (c x 5t x + C y 5t y ) , (6) 

where t' is the corresponding WMAP CMB temperature map and the coefficients c x and c y 
can be determined by minimizing t 2 dean . Using the standard IDL program "regress" we get 
c x = —0.35, c y = —0.78 for the WMAP7 year-1 Ql-band map. We also model and remove 
the pseudo-dipole signal from the released WMAP7 year-1 to year-7 maps of Ql, Q2, VI, 
V2, Wl, W2, W3 and W4 bands, separately^. From the clean maps we calculate their power 
spectra and residual quadruples. Table 1 lists the obtained residual quadrupoles of different 
bands. The overall average clean quadrupole power for all bands is found to be ~ 17.1/xK , 
only ~ 14% of what released by the WMAP teano, indicating that most of the published 



3 In calculation, after adding a deviation to the LOS vector, we always re-scale the LOS vector to keep it 

unitary. 

4 Our map-making codes are publicly released on the website of Tsinghua Center for Astrophysics at 



http://dpc.aire.org.cn/data/wmap/09072731/rclease_vl/source_code/vl/ and on the CosmoCoffce Forum 



at http://cosmocoffcc.info/viewtopic. php?p=4525#4525. 

5 For example, if the two templates 5t x and 6t y are derived from WMAP7 ycar-1 Ql 
band, then they are used to fit the WMAP7 year-1 Ql band single year CMB tem- 
perature map. The WMAP7 single year CMB temperature maps can be found at 



http://lambda.gsfc.nasa.gov/product/map/dr4/maps_forered-da-r9_i_lyr_get.cffh] 



6 The WMAP pseudo-Cl power spectrum for the quadrupole can be found in the WMAP1 release (123 /jK 2 
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Table 1: The residual quadrupole power (/uK 2 ) 



Wave Band 


Ql 


Q2 


VI 


V2 


Wl 


W2 


W3 


W4 


Residual Quadrupole 


16.0 


13.6 


10.4 


13.7 


27.5 


20.2 


12.6 


23.3 


Standard Deviation 


7.0 


3.6 


3.9 


2.8 


7.0 


8.5 


5.3 


27.4 



Note. — For each band, the average residual quadrupole (line 2) and the sample standard deviation (line 
3) is calculated from the 7 single year results. 

WMAP CMB quadrupole is artificial. 



4. DISCUSSION 

We have found that at least three systematical effects inducing pseudo-dipole signal in 
WMAP raw differential temperature data can be summarized into a uniform description: 
the equivalent differential direction error An. There could be more sources with similar 
attributes to these three, but they can all be covered equivalently by a simple An. For 
example, the LOS is determined by observing the Jupiter; however, for each observation, 
the relative direction of Jupiter provides no more than three conditions: (x,y,z), but we 
need four conditions to definitely determine the attitude of the spacecraft: (x, y, z) and rota- 
tion. What's more, the LOS is determined together with the antenna main beam response; 
therefore, the uncertainty of the beam response is also mixed into the LOS uncertainty. 
Moreover, even the WMAP team themselves have found that some systematical effect can 
produce quadrupole- like d eviations on the final CMB map (Jarosik et al. 2007, Fig. 3). 



Recently, iRoukemal (120101 ) pointed out that a small timing error during the step of calibra- 
tion of the raw TOD could correspond to adding a pseudo-dipole difference signal with no 
effect on positional data. It is impossible to exactly estimate all such effects, especially when 
some of them could be unknown for now, but it is easy to model the impact of An upon the 
final CMB map and remove it by model fitting. This is apparently much more feasible and 
effective than estimating and removing all possible sources one by one. 

What is important for model fitting is the feature of a template map, not the absolute 
amplitude. For example, the WMAP foreground removal will be the same even if the three 



for quadrupole) at http://lambda.gsfc.nasa.gov/data/map/powspec/wmap_binned_tt_powspec_yrl_vlpl.txt. 
Although the WMAP team never gave a pseudo-Cl quadrupole in later releases, we have tested and made 
sure that this value is almost the same in all releases from WMAP1 to WMAP 7. 



-7- 



foreground templates ta, t H and t s are all doubled. In this work, the result of removing 
pseudo-dipole signal is not dependent to what amplitude is assumed for the direction error: 
larger amplitude leads to lower fitted coefficient in model fitting and the final correction to 
the CMB map is the same. Here it is worth to mention that the characteristic feature of 
template map of the ps eudo-dipole signal indu ced by a pointing or timing error obtained 
from our previous work (ILiu. Xiong fc Lil 120101) or shown in Fig. 1 of this work has alread y 
been confirmed by other independent works (IMoss. Scott &: Sigurdsonll20ld ; lRoukemall2010l ). 



Our result demonstrates that a pseudo-dipole field can produce deviations in the sky 
temperature map with a structure very similar to the CMB quadrupole pattern published 
by the WMAP team. The template map of pseudo-dipole induced temperature deviation 
shown in Fig. 1 are generated from a differential dipole field which is completely deter- 
mined by the spacecraft velocity and overall-equivalent error of differential direction without 
using any CMB signal. The pseudo-dipole signal in the differential datum for a WMAP 
observation distorts the temperatures for corresponding sky pixels varying observation by 
observation. Therefore, the pseudo-dipole induced temperature map is generated by the 
pseudo-dipole signal combined with the WMAP sc an pattern, which is highly relative to the 
ecliptic plane fainshaw et~ai1l2007l : Ilj et aljlioooj ). The WMAP CMB quadrupole compo- 
nent being highly aligned and close to the ecliptic plane is a long time puzzle in cosmology, 
which now can be naturally explained by the pseudo-dipole signal effect. Keeping insist that 
the published WMAP quadrupole is really cosmological origin now becomes more difficult: 
it is needed the primordial density fluctuations not only occasionally being laid down in the 
plane of the solar system, but also occasionally having almost the same phase as the scan 
strategy of the WMAP mission! 
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